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Abstract 

The  long  term  objective  of  this  program  has  been  to  develop  the  requisite  technologies  for 
integrating  different  sensors  and  actuators  in  desired  combinations  with  electronics  and 
telemetry  on  the  same  substrate.  This  has  been  a  multi-year  program.  We  are  now 
completing  the  program 

Our  concentration  has  been  on  wall-shear-stress  sensors.  We  have  produced  and  calibrated 
floating  element  wall-shear-stress  sensors,  both  active  and  passive.  We  have  developed  a 
calibration  channel  and  a  calibration  procedure  based  on  continuum  isothermal 
compressible  channel  flow.  We  have  also  developed  and  tested  a  wall-shear-stress  sensors 
with  integrated  electronics  initially  based  on  the  Analog  Devices  Inc.  (ADI)  airbag 
accelerometer  circuitry,  and  more  recently  on  one  developed  at  CWRU. 

This  most  recent  work  continued  the  emphasis  on  the  accuracy  of  the  calibration  process 
toward  reducing  the  uncertainty  of  measurement  from  the  past  estimate  of  5%  to  a  goal  of 
1%.  Our  results  are  reported  in  a  M.S  Thesis  recently  completed,  followed  by  some 
additional  Navier-Stokes  computations.  A  journal  manuscript  is  in  preparation. 

Past  Work 

Prior  to  the  start  of  this  grant,  we  built  a  new  calibration  channel  satisfying  the 
requirements  for  non-leakage,  non-choking,  and  having  appropriate  pressure  and 
temperature  measurements.  The  pressure  measurements  are  required  to  evaluate  the 
shear  stress  and  the  temperature  measurements  are  to  assure  the  isothermal  character  of 
the  flow.  It  has  been  shown  in  our  past  work  that  in  order  to  separate  shear  forces  on  the 
sensor  element  from  pressure  forces  on  the  edge  of  the  element,  it  is  necessary  to  know 
the  pressure  gradient  in  the  flow  at  the  sensor  location  both  during  calibration  and  in 
operation. 


Recent  Work 

1.  Improving  accuracy  of  the  calibration  process 

In  an  effort  to  improve  the  accuracy  and  reduce  the  uncertainty  in  the  calibration  process, 
each  element  of  the  procedure  was  investigated  in  detail.  The  basic  calibration 
relationship  is 


r  =  -  (b/2)(dp/cbc)(l-  yM2) 

where  b  is  the  height  of  the  channel,  dp/dx  is  the  measured  pressure  gradient  at  the  sensor 
location,  y  is  the  ratio  of  specific  heats  of  the  test  gas  (nitrogen)  and  M  is  the  Mach 
number  at  the  sensor  location.  We  have  had  our  rotameters  recalibrated  by  the 
manufacturer  so  that  our  Mach  number  determination  is  within  1%.  The  pressure 
measuring  circuits  have  been  cleaned  out  and  now  give  reliable  readings.  Since  it  is  the 
pressure  gradient  that  is  required,  the  absolute  value  of  pressure  is  of  lesser  importance, 
but  we  are  trying  to  get  an  independent  verification  of  the  pressure  transducer  calibration. 


The  largest  inaccuracy  is  in  the  measurement  of  the  channel  height,  b.  Since  this  height  is 
if  the  order  of  200  pm,  it  is  not  measurable  within  1%  by  any  standard  means.  It  is  difficult 
to  bring  a  microscope  to  the  channel  or  vice-versa. 

We  have  instead  used  the  microwave  technique  called  “Microstrip”  (Pozar,  David  M.: 
Microwave  Engineering,  2nd  Ed.,  Wiley,  1998,  pp.160  ft).  A  gold  microstrip  line  is 
printed  on  the  glass  cover  plate  of  the  channel.  The  metal  bottom  of  the  channel  is 
grounded.  The  transmission  frequency  that  minimizes  the  reflected  signal  is  directly 
related  to  the  gap  between  the  glass  and  the  metal.  Our  measurements  indicate  a 
sensitivity  of  about  ±  3pm  which  is  about  ±1 .5%  for  a  channel  height  of  200  pm.  We 
believe  our  overall  uncertainty  now  to  be  about  ±  2%. 

These  results  are  presented  in  the  M.S.  Thesis  of  Mr.  Mehul  Patel  entitled  “Development 
of  MEMS  Wall  Shear  Stress  Calibration  System.”  A  copy  of  the  draft  of  a  journal  article 
based  on  this  work  is  enclosed. 

2.  Navier-Stokes  study  of  flow  around  the  sensor  element 

The  calibration  relations  used  to  date  assume  that  the  flow  about  the  sensor  element  is 
subject  to  the  same  pressure  gradient  as  in  the  center  of  the  calibration  channel.  These 
assumptions  were  examined  this  past  year  in  an  unsponsored  work  by  a  M.S.  student,  Ms. 
Yng-Ru  Chen.  A  presentation  of  this  work  is  planned  for  the  APS/DFD  Annual  Meeting 
in  November  2001 . 

A  two-dimensional  commercially-available  two-dimensional,  incompressible  Navier- 
Stokes  code  was  applied  to  this  problem.  The  flow  domain  included  the  full  channel  height 
plus  the  flow  around  and  beneath  the  sensor  element  and  extended  from  about  one 
channel-height  upstream  -  to  about  a  channel-height  downstream  of  the  sensor  element. 
The  incoming  flow  in  the  main  channel  was  assumed  to  be  fully  developed  laminar  channel 
flow  at  a  Reynolds  number  of  600.  Extra  grid  was  concentrated  in  all  the  passages  around 
and  below  the  sensor  element. 

The  results  show  that  the  pressure  gradient  on  the  channel  side  of  the  sensor  element  was 
essentially  that  in  the  free  stream.  Below  the  sensor  element,  the  pressure  gradient  was  less 
than  in  the  free-stream  by  the  amount  of  the  pressure  losses  in  the  entrance  and  exit 
passages  to  the  channel  underneath  the  element.  The  flow  beneath  the  sensor  element 
displayed  a  parabolic  profile  of  the  appropriate  speed  to  the  pressure  gradient  and  channel 
height.  Because  the  channel  height  beneath  the  sensor  element  was  0.01  time  the  basic 
channel  height,  the  Reynolds  number  beneath  the  sensor  element  is  about  0.0006.  The  fron 
and  rear  edge  surfaces  of  the  sensor  element  displayed  very  unusual  pressure  distributions 
corresponding  to  the  very  low  Reynolds  number  flow  about  them.  The  average  pressure 
on  the  upstream  edge  was  higher  than  in  the  free  stream  immediately  above  that  edge, 
while  the  average  pressure  on  the  downstream  edge  was  correspondingly  lower  than  in  the 
immediately  above  free  stream.  An  analysis  of  these  results  shows  that  use  of  the  assumed 
calibration  relation  systematically  overestimates  the  surface  shear  stress  by  about  0.9%. 
This  factor  has  to  be  introduced  into  the  calibration  process. 


Also  to  be  noted  is  that  the  flow  in  the  passages  around  and  below  the  sensor  element  are 
symmetric,  as  expected  for  this  very  low  Reynolds  number  flow. 
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Development  of  MEMS  Wall  Shear  Stress  Calibration  System 


Mehul  P.  Patel*  and  Eli  Reshotko' 

Case  Western  Reserve  University ,  Cleveland,  Ohio  441-06 

A  new  calibration  procedure  for  the  development  of  wall  shear  stress  sen¬ 
sors  based  on  continuum  isothermal  compressible  flow  has  been  developed. 
The  calibration  relation  is  derived  from  the  momentum  theorem  of  fluid  me¬ 
chanics.  Microelectromechanical  based  floating  element  sensors  for  direct 
measurement  of  wall  shear  stress  are  chosen  for  calibration.  The  calibration 
equation  states  a  linear  relation  between  the  channel  height  and  the  cali¬ 
brated  shear  stress,  resulting  in  an  overall  calibration  uncertainty  derivative 
of  measurement  inaccuracy  from  the  channel  height.  A  novel  microstrip 
technique  formulated  from  the  basic  principles  of  microwave  engineering  has 
been  applied  to  measure  the  channel  height  with  an  increased  level  of  accu¬ 
racy.  The  uncertainties  involved  in  the  measurement  of  all  the  parameters 
that  are  related  in  any  way  to  the  calibration  equation  are  identified,  meas¬ 
ured,  weighed  and  integrated  into  the  equation  to  obtain  a  complete  error 
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analysis  of  the  calibration  procedure.  The  calibration  of  shear  stress  sensors 
is  successfully  completed  with  an  overall  accuracy  of  2%. 


I.  Introduction 

The  study  of  drag  caused  by  a  fluid  flow 
over  bodies  of  varied  shapes  is  of  vital  im¬ 
portance  to  aerodynamic  and  hydrodynamic 
design.  Knowledge  of  skin  friction  is  of 
paramount  importance,  because  friction  de¬ 
termines  a  major  part  of  the  drag  and  can  be 
an  important  factor  in  limiting  or  affecting 
the  performance  of  a  vehicle  or  structure. 
Wall  shear  stress  is  a  quantity  that  is  gener¬ 
ally  inferred  since  it  cannot  be  easily  meas¬ 
ured.  Despite  the  considerable  work  done  on 
measurement  of  wall  shear  stress,  there  is 
still  a  need  for  an  easily  implementable  and 
reliable  technique  of  direct  measurement  of 
wall  shear  stress.1,  2'  3'  4  The  indirect  tech¬ 
niques  work  under  assumptions  that  relate 
the  shear  stress  to  the  quantities  that  are  di¬ 
rectly  measured.  The  only  direct  measure¬ 
ment  technique  is  the  floating-element  tech¬ 


nique,  which  use  floating  elements  that  dis¬ 
place  against  spring  suspension  systems 
when  exposed  to  shear  stress.  Previous  work 
describes  the  design,  fabrication,  testing  and 
calibration  of  shear  stress  sensors  based  on 
the  floating-element  method  with  an  overall 
inaccuracy  of  5%.5,6, 7' 8  A  flow  channel  that 
could  sustain  a  two-dimensional  continuum 
flow  of  nitrogen  gas  was  designed  for  the 
calibration  of  Microelectromechancal  sys¬ 
tems  (MEMS)  sensors.  An  expected  linear 
relationship  between  sensor  deflection  and 
shear  stress  was  demonstrated  and  correlated 
with  mechanical  modeling. 

This  paper  presents  the  refinement  in  the 
development  of  wall  shear  stress  calibration 
system  based  on  continuum  isothermal 
compressible  flow  theory.9  Flow  channel  is 
calibrated  prior  to  the  actual  calibration  of 
MEMS  shear  stress  sensors.  A  novel  micro- 
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strip  technique  based  on  the  principles  of 
microwave  engineering  is  formulated  to  re¬ 
duce  the  largest  sources  of  error  identified 
from  channel  height  measurements.  Finally, 
a  complete  error  analysis  of  the  calibration 
procedure  has  been  carried  out  and  inte¬ 
grated  into  the  calibration  equation. 

II.  MEMS  Wall  Shear  Stress  Sensors 

The  sensors  used  for  calibration  in  the 
presented  work  are  designed  and  fabricated 
using  MEMS  techniques.  The  floating  ele¬ 
ment  design  on  the  MEMS  shear  stress  sen¬ 
sors  are  based  on  a  lateral  resonant  structure, 
a  design  that  has  been  used  in  a  number  of 
microsensing  and  microactuating  devices.10, 

12  Lateral  resonant  structure  features  and 
theory  are  detailed  in  Ref.  7.  A  square  float¬ 
ing  element  with  tether  suspension  beams 
are  microfabricated  in  polysilicon.  Located 
at  the  center  of  the  device,  the  floating  ele¬ 
ment  is  suspended  2.5  microns  above  the 
lower  electrodes  by  the  folded  beam  suspen¬ 
sion.  The  deflection  of  the  folded  suspension 


beam  with  applied  shear  force  is  linear.  A 
picture  of  a  MEMS  sensor  developed  at 
Case  Western  Reserve  University  (CWRU) 
and  Advanced  MicroMachines  for  Industry 
(AMMI)  used  for  calibration  is  shown  in 
Fig.  1.  The  sensors  use  floating  elements 
that  displace  against  spring  suspension  sys- 


Fig.  1  CWRU/AMMI  floating  element  shear 
stress  sensor. 


terns  when  exposed  to  shear  forces,  thus 
qualifying  as  a  direct  measurement  tech¬ 
nique. 

The  application  of  MEMS  technologies 
to  fluid  mechanics  offers  new  possibilities 
for  improved  understanding  of  the  physics 
of  wall  boundary  layers.13,  14,  15  In  order  to 
perform  a  complete  analysis  of  wall  shear 
stress,  fluctuations  occurring  even  in  the  mi- 


cro  scale  (small  eddies)  have  to  be  resolved. 
The  scales  of  motion  in  a  turbulent  boundary 
layer  are  of  the  order  of  few  hundred  mi¬ 
crons  and  the  scale  of  the  smallest  eddies  is 
about  100pm  while  the  lifetime  of  these 
structures  of  the  turbulent  boundary  layers  is 
in  milliseconds.16  The  sensor  dimensions 
thus,  needs  to  be  smaller  than  the  size  of  the 
smallest  eddies.  The  most  important  features 
for  a  sensor  that  is  used  for  measuring  small- 
scale  shear  stress  fluctuations,  are  small  ac¬ 
tive  sensor  area  and  fast  response  time.  Us¬ 
ing  micro-machined  shear  stress  sensors, 
high  spatial  and  temporal  resolution  can  be 
obtained  for  the  measurement  of  wall  shear 
stress.  Arrays  of  such  sensors  would  allow 
the  accurate  mapping  and  analysis  of  turbu¬ 
lent  structures,  with  speeds  sufficient  to  de¬ 
tail  their  formation  and  destruction. 

III.  Theoretical  Analysis 

The  schematic  of  the  flow  channel, 
which  is  used  for  sensor  calibration,  is 
shown  in  Fig.  2.  The  flow  inside  the  high 


aspect  ratio  (96-145)  continuum  flow  chan¬ 
nel  is  approximated  to  be  a  two-dimensional 
flow.  A  steady  continuous  flow  of  nitrogen 
gas  is  supplied  from  a  regulated  nitrogen 
cylinder.  For  the  channel  with  uniform  cross 
sectional  area,  it  is  assumed  that  the  pressure 
variation  is  only  along  the  flow  direction. 
The  momentum  theorem  of  fluid  mechanics 
is  used  to  derive  the  relation  that  is  used  for 
the  calibration  of  the  shear  stress  sensors. 


Fig.  2  Schematic  of  the  flow  channel  base,  shim  and  fi¬ 
ring. 

which  makes  the  relation  independent  of 
whether  the  flow  is  laminar,  transitional  or 
turbulent. 

Assumptions  for  the  fluid  flow  are: 

1 .  The  flow  is  steady. 
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2.  The  local  velocity  or  Mach  number  dis¬ 
tributions  can  be  represented  by  their 
mean  values  across  the  channel  height. 

3.  Temperature  is  constant. 

A  control  volume  defined  as  an  imagi¬ 
nary  fixed  volume  through  which,  a  fluid 
may  flow  is  shown  in  Fig.  3.  The  pressure 
gradient  inside  the  flow  channel  is  the  driv- 

O 

ing  force  of  the  fluid. 
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Fig.  3  Control  volume  of  the  flow  channel. 

A.  Continuity  Equation 


The  continuity  equation  states  that 


—  AmJyRT  =  constant  (2) 
RT 


In  differential  form  (2)  becomes. 


1  dT  .  dp  .  dA  1  dM2 

- 1 - 1 - 1 - - 

2  T  p  A  2  M- 


=  0 


Since  cross-sectional  area  is  constant  and 
gas  is  assumed  to  be  isothermal,  i.e. 


dA 

A 


dT_ 

T 


=  0, 


Therefore,  dM  2  =  -2  M 2 


dp 


P 


(3) 


B.  Momentum  Balance 

Applying  momentum  balance  on  a  rigid 
control  volume  of  differential  size,  we  get: 


pAV  -  constant  (1) 


\pV 2  xbx w]1+(iv  -[ pV2  x&xw],  = 


For  an  ideal  gas,  p  =  pRT  and  V  =  Ma , 

[pxbxw  -  (p  +  dp)x  b  x  w]-  2r  Wwdx 

where,  a  =  -JyRT  ,  Therefore  (1)  implies, 


5 


^E.+—(pv2) 

dx  dxy  ' 


dp  d 
—  +  — 
dx  dx 


-?-M2yRT 

RT 


b 
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dp  d 
—  +  — 
dx  dx 


(pyM-) 


b 

T"=-9 


dp  dp  dM 2 

-J-  +  yM-  —  +  yp—- 
dx  dx  dx 


(4) 


From  (3), 


dM 2  2  dp 

yp——  =  -2yM  -f 
dx  dx 


(5) 


Substituting  (5)  in  (4), 


b 

T"'  “  2 


^(i -*') 

dx 


(6) 


where,  b  is  the  channel  height,  TM.is  the  wall 


.dp  . 

shear  stress  and  —  is  the 
dx 


pressure  gradient 


at  the  sensor  location.  The  Mach  number 
correction  term  in  (6)  incorporates  the  com¬ 
pressibility  effect  for  an  isothermal  channel 


Phigh 


Piow- 


Fig.  4  Effective  stresses  on  floating  element. 


flow  in  the  calibration  equation. 

There  are  other  forces  acting  on  the 
floating  element,  which  are  also  considered 
for  calculating  the  shear  stress  on  the  ele¬ 
ment.  Fig.  4  shows  the  details  of  the  floating 
element  within  the  bed.  The  pressure  gradi¬ 
ent  generated  inside  the  flow  channel  is  the 
driving  force  of  the  fluid.  A  “gauge  thick¬ 
ness  factor”  as  shown  in  Fig.  4  arises  due  to 
the  pressure  differential  across  the  element, 
which  forces  the  element  to  deflect  in  the 
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direction  of  the  lower  pressure.  This  addi¬ 
tional  gauge  force  is  acting  on  the  floating 
element  and  needs  to  be  subtracted  in  order 
to  accurately  determine  the  wall  shear  stress. 

F gauge,  can  be  measured  as  a  pressure  dif¬ 
ferential,  A P gauge,  on  the  fore  and  the  aft 
edge  surfaces  of  the  element.  Let  the  dimen¬ 
sions  t  and  a,  be  the  thickness  and  the  width 
normal  to  the  page,  respectively.  The  gauge 
factor  F smiSe,  can  then  be  written  as, 


F gauge  =  A  P gauge 


Let  V  be  the  length  of  the  floating  ele¬ 
ment.  The  gauge  pressure  difference  A Pga„ge, 
is  the  pressure  differential  acting  across  the 
length  c  of  the  element  and  is  equal  to 
dp h 


-c- 


dx 


.  Therefore, 


F gauge  (ta)(  C  dy,  ) 


the  floating  element  by  the  top  surface  area 
of  the  element  ‘ac': 


x 


gauge 


(7) 


Additionally,  there  is  a  small  quantity  of 
fluid  flow  underneath  the  floating  element 
that  contributes  to  yet  another  shear  stress 
on  the  lower  surface  of  the  floating  element, 
which  is  comparable  to  the  gauge  stress  but 
very  small  compared  to  the  top  surface  shear 
stress.  Let  ‘g’  be  the  gap  spacing  between 
the  floating  element  and  the  fixed  substrate. 
Then  this  shear  stress  is  expressed  as, 


g_dp_ 
2  dx 


(8) 


The  effective  applied  shear  stress  after 
completing  the  force  balance  equation  in  the 
steady  state  can  be  written  as, 


The  effective  gauge  stress  Tgauge,  is  then 
calculated  by  dividing  the  forces  acting  on 


r  „  =  t 

eJ7  v 


+  r  +r 

gauge  gap 


1 


(6),  (7)  and  (8)  give: 


Xeff 


dp_ 

dx 


-(l-yM2)  +  ^-  +  f 


(9) 


The  above  equation  shows  the  total  ap¬ 
plied  force  acting  on  the  floating  element 
divided  by  the  upper  surface  area  of  the 
element  during  the  actual  calibration  of  the 
sensor.  The  shear  stress  on  the  working  sur¬ 
face  of  the  floating  element,  rw,  is  given  by 


T.  t  eff  ^  pane  ^  pap  ^  eff 


dp 

~dx 


7  g  ^ 

—  +  t 
2  j 

(10) 


IV.  Calibration  Apparatus 

The  use  of  a  floating  element  sensor  to 
measure  rK.  requires  a  simultaneous  meas¬ 
urement  of  the  local  pressure  gradient.  Sen¬ 
sor  calibration  involves  the  critical  job  of 
relating  the  sensor  output  to  a  primary  shear 
stress  standard.  The  two-dimensional  iso¬ 


thermal  compressible  constant-area  channel 
flow  is  regarded  as  such  a  standard. 

A.  The  Flow  Channel 

A  three-dimensional  schematic  of  the  as¬ 
sembled  flow  channel  is  shown  in  Fig.  2. 
The  channel  consists  of  a  base  plate  made  of 
stainless  steel,  a  glass  cover  plate  and  a  plas¬ 
tic  shim.  The  thin  plastic  shim  that  defines 
the  height  and  width  of  the  channel  is  placed 
on  the  top  of  the  base  plate.  The  assembled 
flow  channel  is  formed  by  clamping  the 
glass  cover  plate  over  the  shim  with  the  help 
of  another  steel  plate  bolted  with  the  bottom. 
The  channel  height  can  be  varied  by  using 
shims  of  different  thicknesses  or  by  tighten¬ 
ing  and  loosening  the  top  screws.  The  chan¬ 
nel  height  is  varied  to  provide  different 
ranges  of  wall  shear  stress  at  the  sensor  loca¬ 
tion  inside  the  channel  for  the  given  flow 
rates.  A  U-ring  is  mounted  outside  the  shim 
to  provide  a  seal  against  side  leakage.  The 
assumption  of  a  two-dimensional  continuum 
fluid  flow  holds  good  due  to  the  high  aspect 


8 


ratio  of  the  channel.  Additional  detail  on  the 
flow  channel  design  is  discussed  in  Ref.  8. 

The  micro-fabricated  shear  stress  sensor, 
mounted  on  a  microchip,  is  placed  at  a  dis¬ 
tance  of  16.94  mm  from  the  channel  exit  and 
is  made  flush  with  the  channel  bed  by  insert¬ 
ing  it  inside  a  machined  groove.  The  groove 
dimension  is  such  that  it  enables  the  top  sur¬ 
face  of  the  chip  to  fall  in  the  same  plane  as 
the  channel  bed. 

There  are  leveling  screws  as  seen  in  Fig. 
5  underneath  the  mounting  plate  of  the  sen¬ 
sor  chip  for  very  fine  adjustments  to  make 
the  sensor  surface  flush  with  the  channel 
bed.  A  groove  with  a  depth  of  0.068  inches, 
which  is  standard  for  l/8lh  inch  O-ring  static 
seals,  is  machined  into  the  base  to  accom- 


Fluid  Flow 


Flow  Channel  Base 


Channel  Exit 


Sensor  Die  Electrical  Lead 


3/56  Fine  Thread 


Fig.  5  Schematic  of  sensor  mounting  apparatus. 


modate  a  U-ring  for  separating  the  glass 
from  metal  and  obtaining  a  leakage  free 


Flow  Channel  Base  (Aluminum) 


- . -Gold  Lines  MM}.  Shim  (Plastic)  tiSS82®ZSiCmtact-PaJ 

Fig.  6  Front  view  of  the  channel 
flow  channel.  Fig.  6  shows  the  location  of 
the  U-ring  in  preventing  the  leakage  from 
the  channel. 

B.  Pressure  Transducer 

There  are  five  miniature  pressure  taps 
that  are  drilled  inside  the  channel.  One  end 
of  each  tap  opens  inside  the  channel  and  the 
other  end  of  the  taps  are  attached  to  different 
hoses,  which  are  connected  to  a  NIST  cali¬ 
brated  pressure  transducer  OMEGA® 
PX236.  The  first  pressure  tap  is  located  at  a 
distance  of  38.735mm  from  the  channel  end 
and  the  fifth  tap  is  at  a  distance  of  8.25mm 
from  the  channel  end.  A  fixed  distance  of 
7.625mm  is  maintained  between  each  of  the 
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taps  to  ensure  uniformity.  These  taps  enable 
direct  measurement  of  the  static  pressures  at 
these  locations.  A  second  order  polynomial 
curve  fit  is  made  once  the  pressure  readings 
at  these  five  known  locations  are  taken.  The 
pressure  gradient  is  obtained  from  differen¬ 
tiating  this  curve  fit.  The  pressure  at  the  sen¬ 
sor  location  is  obtained  by  interpolation  in 
the  pressure  distribution. 

C.  THERMOCOUPLES 

Continuous  monitoring  of  the  tempera¬ 
ture  at  different  locations  inside  the  flow 
channel  is  necessary  to  ensure  isothermal 
conditions  inside  the  channel.  For  this  rea¬ 
son  three  tungsten  wires  acting  as  P-type 
thermocouples  have  been  inserted  in  1/1 6th 
of  an  inch  diameter  holes  at  three  different 
locations  inside  the  channel  bed  from  the 
bottom  through  drilled  holes.  The  locations 
of  these  probes  are  at  distances  of  8.25  mm, 
23.49  mm  and  38.735  mm  respectively  from 
the  channel  end.  The  thermocouple  tips  are 
exposed  to  the  inside  of  the  channel.  The 


thermocouple  probes  are  inserted  in  such  a 
manner  that  they  provide  a  real  time  tem¬ 
perature  measurement  at  all  three  locations 
without  causing  any  obstruction  to  the  fluid 
flow. 

D.  Rotameter 

It  is  necessary  to  maintain  a  continuous 
non-fluctuating  fluid  flow  for  the  calibration 
of  the  sensor.  The  rotameter  controls  and 
measures  the  flow  rate  that  is  required  to 
calculate  the  Mach  number  of  the  fluid  flow, 
which  is  eventually  required  for  the  calibra¬ 
tion  of  the  shear  stress  sensors.  Two  rotame¬ 
ters  of  different  ranges  of  flow  rates,  certi¬ 
fied  with  NIST  calibration  standards,  are 
used  in  the  experiment. 

1 .  Low  flowrate:  288.0-7590.0  (std  ml/min) 

2.  High  flowrate:  1384.6-46132.7  (std 
ml/min) 

Both  rotameters  are  certified  with  NIST 
calibration  to  ensure  their  accuracy. 
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V.  Floating  Element  Shear  Stress 
Micro-sensor  with  Differential  Ca¬ 
pacitive  Read-out  Scheme 

A  differential  capacitive  read-out 
scheme  has  been  adopted  to  measure  the  po¬ 
sition  change  of  the  floating  element.  When 
subjected  to  fluid  flow,  the  shear  force  dis¬ 
places  the  floating  element,  causing  the  ele¬ 
ment  electrodes  to  move  closer  to  one  set  of 
fixed  electrodes  and  away  from  another.  The 
two  capacitances  formed  by  the  movable 
and  fixed  electrodes  change  accordingly;  C] 
increasing  and  C?  decreasing  are  then  con¬ 
verted  to  voltages  and  amplified.  On  chip 
circuitry  on  the  sensor  chin  provides  force- 
balance  sensor  operation.17,  18  Floating  ele¬ 
ment  voltage  signal  is  modulated  by  an  ex¬ 
ternal  referenced  oscillator  and  fed  into  a 
preamplifier.  This  signal  is  fed  back  to  the 
stationary  electrodes  as  a  bias  signal.  This 
bias  incites  an  electrostatic  force  to  attract 
the  floating  element  towards  the  far  elec¬ 
trode  and  away  from  the  near  electrode. 


Electrostatic  forces  in  these  devices  are  at¬ 
tractive  only,  such  that  the  net  attractive 
force  in  a  given  direction  of  interest  causes 
the  floating  element  to  move  back  to  its  zero 
displacement  position.  Details  of  the  differ¬ 
ential  capacitive  read-out  scheme  for  MEMS 
shear  stress  sensors  used  in  this  work  are 
presented  in  Ref.  8. 

VI.  Measurement  of  Parameters 

The  parameters  that  need  to  be  measured 
experimentally  for  shear  stress  determina¬ 
tion  from  the  calibration  relation  (9)  are: 

1.  Channel  height  ‘b\  as  defined  by  the 
thin  shim  thickness 

2.  Pressure  gradient  ,  at  the  sensor 
location 

3.  Mach  number  ‘M\  at  the  sensor  location 

4.  Thickness  V’  of  the  floating  element  and 

5.  The  gap  ‘ g ’  between  the  floating  element 
and  the  electrostatic  shield. 
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A.  Channel  Height 

In  the  series  of  experiments  that  were 
performed  for  the  calibration,  the  channel 
heights  ranged  from  180|im  to  260p.m.  The 
height  is  varied  depending  on  the  dynamic 
range  of  the  shear  stress  sensors  that  are 
used  for  calibration.  Since,  a  linear  depend¬ 
ency  between  channel  height  and  shear 
stress  exists  in  the  calibration  relation  as 
shown  in  (9),  the  percent  error  in  the  meas¬ 
urement  of  channel  height  will  contribute 
the  same  percentage  error  in  the  calculation 
of  shear  stress.  Consequently,  channel 
height  b,  the  most  difficult  parameter  to 
measure  amongst  all  five  parameters,  is  also 
the  most  important  one.  In  the  previous 
calibration  technique,  the  channel  height 
was  measured  using  an  optical  microscope. 
The  uncertainty  in  the  measurement  of  the 
height  with  the  optical  measurement  tech¬ 
nique  was  estimated  to  be  ±20pm.  In  the 
present  work,  a  novel  micro-strip  technique 
has  been  designed  from  the  fundamental 


principles  of  microwave  engineering  and  has 
been  demonstrated  to  reduce  uncertainty  in 
the  measurement  of  channel  height. 

Al.  Micro-Strips 

Micro-strip  lines  are  thin  planar  trans¬ 
mission  lines,  which  support  electromag¬ 
netic  waves  at  microwave  frequencies. 
These  lines  are  fabricated  by  photolitho¬ 
graphic  processes,  in  which  conductive  ma¬ 
terials  like  gold  or  copper  are  deposited  on  a 
substrate.  The  geometry  of  such  a  micro¬ 
strip  line  is  shown  below.  A  conductor  of 
width  w  is  printed  on  a  thin,  grounded  di¬ 
electric  substrate  of  thickness  d  and  relative 
permittivity  er .  A  sketch  of  the  field  lines  is 
shown  in  Fig.  7. 


Fig,  7  Microstrip  transmission  line,  (a)  Geome¬ 
try.  (b)  Electric  and  Magnetic  field  lines. 
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Using  a  quasi-static  approach,  the  phase  The  ‘ w/d ’  ratio  is  given  by 

velocity  v  and  propagation  constant  /3  can 

be  approximated  as,  33L-Z. 

b  ft 
c 

Phase  velocity,  vp  =  —j= 

Propagation  constant,  (3  -  koy[E^ , 

for  w/b  >2  (12) 


B-l-  ln(2fi  -l)+ 


where,  c  =  3xl08  m/sec  is  the  speed  of 
light  and  Ee  is  the  effective  dielectric  con¬ 
stant  of  the  micro-strip  line.  For  a  lossless 
line,  the  characteristic  impedance  Z0,  which 
is  the  ratio  of  voltage  to  current  for  the  trav¬ 
eling  wave,  is  expressed  as,Z0  = 

where  L  and  C  are  inductance  and  capaci¬ 
tance  per  unit  length  of  the  line.  Given  the 
substrate  (air)  thickness,  b,  and  conductor 
width,  vv,  the  characteristic  impedance  of  the 
line,  Z0,  is 

_  _ 120ft  _ 

°  ”  + 1 393 + °-66v  ^Wb + 1  -444)- 

for  w/b  >1  (11) 


377ft 

Where’B=^ 


A2.  Design  of  the  Gold  Line 

A  gold  line  of  width  w  and  length  /,  as 
shown  in  Fig.  8,  is  deposited  on  the  lower 
surface  of  the  glass  cover  plate  by  a  shadow 
mask  process.20,  21  The  shim  that  separates 
the  glass  cover  plate  from  the  steel  base 
plate  defines  the  channel  height  b.  The  me¬ 
dium  between  the  plates  is  air.  The  effective 
dielectric  constant  of  air,  Ee ,  is  1.  Now,  if 
the  characteristic  impedance,  Z0 ,  and  width 
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Fig.  8  The  length  of  the  transmission  line. 


w  are  fixed,  the  distance,  b  can  be  calculated 
from  (11)  and  (12). 

In  the  present  design,  an  electromagnetic 
RF  wave  from  a  source  inside  the  Network 
Analyzer-HP4395A  is  sent  down  the  gold 


line.  The  characteristic  impedance  of  the 
transmission  line  is  50£2.  The  purpose  of 
fixing  the  characteristic  impedance  to  50Q  is 
because  of  its  standard  usage  in  all  types  of 
applications.  Fig.  9  is  the  block  diagram  of 
the  set  up  for  the  channel  height  measure¬ 
ment.  Fig.  10  is  the  analogy  diagram  to  sim¬ 
plify  and  deduce  the  fact  that  the  present  de¬ 
sign  is  no  different  than  the  actual  applica¬ 
tion  after  taking  into  effect  the  effective 
permittivity  of  the  glass  and  air. 

(11)  and  (12)  show  that  the  interaction  of 
the  micro-strip  with  the  conducting  plane  is 
extremely  sensitive  to  their  separation  dis- 


Fig.  9  Channel  height  measurement  setup  using  HP4395A  -  Network  Analyzer. 
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Stainless  steel 


glass 


1 


Fig.  10  Microstrip  application  -  channel  height  measurement  analogy  diagram. 


tance  (channel  height).  The  channel  base  is 
treated  as  an  infinite  conducting  plane  some 
distance  away  from  the  transmission  line.  So 
by  measuring  the  reflection  and  the  trans¬ 
mission  effects,  the  distance  could  be  deter¬ 
mined.  The  value  of  effective  dielectric  con¬ 
stant,  ee  which  is  1,  and  the  characteristic 

impedance  Z0  that  is  preset  at  50£2  are  sub¬ 
stituted  in  (11)  and  (12)  keeping  a  fixed 
channel  height  b,  say  200pm.  The  w/b  ratio 
is  4.9  for  the  given  values  of  Z0  and  b.  The 
width  of  the  gold  line  w  obtained  from  the 
w/b  ratio  that  is  equal  to  983pm.  Length  of 
the  gold  line  is  also  an  important  parameter 
to  deduce  since  the  reflection  coefficient  de¬ 
pends  heavily  on  it. 


The  idea  is  to  obtain  the  reflection  coef¬ 
ficient  equal  to  zero,  which  is  possible  only 
in  a  special  case  where  the  total  signal  that  is 
transmitted  is  completely  reflected  by  180° 
out  of  phase.  Thus  when  the  length  is  A/2, 
the  reflected  wave  which  is  completely  out 
of  phase,  cancels  out  the  transmitted  wave 
and  it  seems  that  there  is  nothing  there  even 
though  everything  gets  back.  A  large  value 
of  return  loss  corresponds  to  a  small- 
reflected  signal  just  as  a  large  value  for  in- 


Gold  line  sputtered  on  glass 


Fig.  11  Design  of  the  transmission  line  (mate¬ 
rial:  gold). 
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sertion  loss  corresponds  to  a  small- 
transmitted  signal.  A  marker  is  used  to  iden¬ 
tify  the  frequency  at  which  the  lowest  reflec¬ 
tion  coefficient  is  obtained.  The  frequency 
sweep/magnify  feature  allows  zooming  in  a 
particular  range  of  operating  frequencies  to 
obtain  a  much  closer  look  at  the  values  of 
the  reflection  coefficient. 

As  shown  in  Fig.  11,  the  length  could  be 
anything  that  satisfies  the  relation  /  =  X/2  + 
n,  where  n  is  a  positive  integer.  The  condi¬ 
tion  of  obtaining  T  =  0  is  satisfied  for  all 
values  of  /  as  far  as  it  satisfies  the  given  rela¬ 
tion.  The  length  is  fixed  to  X/2  in  this  case 
because  of  the  restraints  from  the  machine 
which  is  used  for  depositing  gold  on  the 
glass  plate.  The  maximum  length  of  the 
glass  that  the  machine  can  handle  is  11.4  cm 
in  length.  So  the  glass  cover  plate  was  sepa¬ 
rated  into  two  pieces  with  one  piece  of 
length  10.16  cm.  The  exact  design  of  the 
gold  line  is  shown  in  Fig.  8.  The  average 
length  of  the  gold  line  is  obtained  by  taking 


the  mean  of  the  inside  line  and  outside  line. 
This  length  is  321.29mm.  The  operating  fre¬ 
quency  is  obtained  from  the  relation,  c  =  ^.f; 

where  X  =  21.  The  frequency  thus  obtained 

is  466.867  MHz.  Since  the  range  of  the  op¬ 
erating  frequencies  for  the  HP  4395A  is 
lOKHz  ~  500  MHz,  the  maximum  wave¬ 
length  X  is  limited  to  600mm.  This  is  the 

reason  for  the  curved  shape  of  the  gold  line, 
rather  than  a  straight  line,  which  will  have 
reduced  losses.  The  higher  the  operating 
frequency,  the  shorter  the  length  of  the 
transmission  line  and  better  the  reflection 
characteristics. 

Figs.  12  and  13  are  the  screen-shots 
taken  from  the  Network  Analyzer  HP4395A 
display  screen  showing  a  minimum  reflec¬ 
tion  occurring  at  458.5  MHz  and  456.7  MHz 
respectively.  The  reflection  coefficient  of 
the  transmission  line  for  the  given  designed 
parameters  should  be  equal  to  zero  in  an 
ideal  case,  but  there  are  return  losses  for  the 
reasons  mentioned  below. 
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1.  The  gold  line  as  shown  in  Fig.  8  has  got 
three  curves  that  contribute  erratic  re¬ 
flection  of  the  signal  in  unwanted  direc¬ 
tions. 

2.  The  launching  of  the  co-axial  from  the 
HP  4395A  to  the  gold  line  involved  a 
mechanism,  which  made  a  compressed 
fit  between  the  co-axial  and  the  line. 
This  method  of  launching  the  co-axial 
scrapes  away  some  percent  of  the  line  at 
every  attempt  to  make  a  connection.  Fur¬ 
thermore,  at  least  8  to  10  such  attempts 
are  required  for  successful  connections 
to  take  a  complete  set  of  data  for  a  par¬ 
ticular  channel  height. 

3.  Low  interference  of  the  electromagnetic 
waves  is  caused  due  to  the  metalization 
on  the  sensor,  which  is  placed  exactly  on 
the  bottom  side  of  the  gold  line  with  its 
face  towards  the  line  and  is  separated  by 
the  distance  of  the  channel  height. 

4.  The  width  of  the  gold  line  is  not  constant 
throughout  the  entire  length  as  shown  in 


Fig.  8.  Instead,  there  are  variations  in  the 
width  at  quite  a  few  places.  The  width 
was  then  made  constant  by  manual 
scraping  of  the  gold  with  a  sharp  blade. 
Such  varying  width  causes  reflection 
losses. 

These  reasons  attribute  toward  a  non¬ 
zero  value  of  the  reflection  coefficient  and 
the  results  drift  by  extremely  small  values  as 
return  losses.  So  the  channel  height  needs  to 
be  backtracked  from  the  obtained  value  of 


511  .  o?  1*3  1.5  i%<  REF  -M3  -12.747  dB 


Fig.  12  Network  analyzer  plot  showing  minimum 
reflection  at  458.5  MHz. 
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Fig.  13  Network  analyzer  plot  showing  minimum  reflection  at 
456.47  MHz. 


return  losses.  The  reflection  coefficient  is 
obtained  in  dB  (decibels)  from  the  HP 
4395A.  These  values  are  then  used  to  obtain 

the  reflection  coefficient  T.  Ref.  19  details 
the  theory  and  the  equations  that  were  used 
for  obtaining  the  reflection  coefficient  T. 

The  plots  from  the  network  analyzer  are  also 
shown  at  the  end.  The  output  from  the  net¬ 


work  analyzer  is  obtained  as  return 
losses  from  which  the  values  of  T 
are  calculated. 

From  Appendix  A.  the  estimate 
for  the  uncertainty  in  the  measure¬ 
ment  of  channel  height  is  ±1.5%. 

B.  Pressure  Gradient 

The  negative  pressure  gradient 

along  the  flow  path  is  charac- 

dx 

terized  with  the  help  of  a  pressure 
transducer  OMEGA®  PX236  that  is 
connected  to  five  pressure  taps, 
which  are  evenly  located  inside  the 
flow  channel.  The  pressure  transducer  is 
calibrated  according  to  the  NIST  standards 
and  procedures.  A  common  switch  is  used  to 
connect  the  pressure  transducer  to  each  of 
the  five  pressure  taps;  thus  eliminating  the 
need  of  more  pressure  sensors.  Once  the 
pressure  readings  are  obtained,  the  dynamic 
property  of  the  flow  can  be  inferred  from  the 
plot  of  the  pressure  gradient.  A  fifth  order 
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polynomial  curve  fit 

p(x )  =  ax5  +  bx 4  +  cxJ  +  dx~  +  ex  +  /  is 
then  made  to  these  pressure  data  to  obtain 
the  pressure  at  the  sensor  location  x.  The 
pressure  gradient 

p  (x)  =  5 ax*  +  4 bx3  +  3 ex'  +  2dx  +  e 

at  the  sensor  location  V  is  obtained  from 
interpolation  of  the  curve.  The  estimated 
uncertainty  of  pressure  gradient  determina¬ 
tion  for  the  presented  work  is  ±0.2%. 

C.  Mach  number 

Mach  number  is  determined  after  meas¬ 
uring  the  average  fluid  velocity  at  the  sensor 
location  Vsensor.  Vsensor  is  a  function  of  the 
volumetric  flow  rate  ‘ Qsensor  and  the  cross 
sectional  area  V\’  of  the  fluid  flow. 


The  flow  rate  Qsens0r  at  the  sensor  loca¬ 
tion  is  determined  from  the  measurements 
taken  upstream  of  the  channel  by  a  NIST 
calibrated  variable  area  flowmeter.  The  scale 
of  the  rotameter  is  calibrated  on  the  assump¬ 
tion  that  its  output  is  at  atmospheric  pres¬ 
sure,  which  actually  is  subjected  to  back 
pressure.  This  pressure  force  pushes  the  ro¬ 
tameter  float  to  slide  down  the  scale,  result¬ 
ing  in  a  measurement  that  is  lower  than  the 
actual  flow  rate  Qrot.  To  eliminate  this  small 
error  in  the  measurement  of  the  flow  rate,  a 
simple  correction  factor  can  be  determined 
for  differing  back-pressures.22  Assuming 
isothermal  flow,  if  pressure  at  the  exit  of  the 
rotameter,  Prot  is  measured,  for  a  given  flow 
rate  reading,  Qroh  the  flow  rate  at  the  sensor 
location,  Qsensor  is  given  by, 
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where,  P sensor  is  the  pressure  at  the  location 
of  the  sensor  and  Palm  is  atmospheric  pres¬ 
sure.  The  velocity  at  the  sensor  location  is 
then  measured  by  dividing  the  actual  flow 
rate  at  the  sensor  location,  Qsenson  by  the 
cross-sectional  area  under  through  the  fluid 
flows. 


Mach  number  is  then  determined  by 
comparing  the  sensor  velocity  to  the  local 
speed  of  sound  ,  Vstmnd  at  the  measured  tem¬ 
perature. 


The  Mach  number  where  choking  occurs 
for  an  isothermal  flow  of  a  perfect  gas  in  a 
constant  duct  area  in  the  presence  of  wall 
friction  is  given  by  1/Vy.23 


D.  Thickness  and  Gap  Measurements 

The  values  of  the  thickness  of  the 
floating  element  ‘t’  and  the  gap  spacing  be¬ 
tween  the  floating  element  and  the  fixed 
substrate  ‘g’  are  obtained  from  the  actual 
design  specifications  of  the  floating  element. 
The  values  of  ‘g’  and  ‘t’  are  0.0025  mm  and 
0.0020  mm  respectively.  Even  a  10%  varia¬ 
tion  from  design  values  of  V  and  ‘g’  con¬ 
tribute  negligibly  to  the  error  in  meas¬ 
urements. 

Vll-A.  Calibration  Flow  Channel  Test 
Results 

Two  sets  of  tests  were  run.  In  the  first, 
the  channel  was  run  without  a  sensor  in 
place.  In  the  second  set,  a  sensor  was  in¬ 
stalled  and  calibration  run  implemented. 

Initial  experiments  on  the  flow  channel 
were  performed  without  the  shear  stress  sen¬ 
sor  in  place  and  were  aimed  at  ascertaining 
the  flow  channel  performance  and  capabili¬ 
ties.  The  pressure  readings  along  with  the 
Mach  number  determinations  were  com- 
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pared  with  the  theoretical  values  for  a  steady 
one-dimensional  isothermal  flow  with  fric¬ 
tion  in  a  constant-area  duct.  The  well-known 
phenomenon  of  choking,  the  limiting  condi¬ 
tion  for  such  a  flow  in  a  constant-area  duct, 
is  also  observed  with  the  present  calibration 
flow  channel. 

Fig.  14  shows  the  nearly  variation  of 
Mach  number  with  increase  in  the  flow  rate 
inside  the  flow  channel  at  channel  height 
230  pm.  Mach  0.40  was  obtained  at  a  flow 
rate  of  43960  seem  at  a  channel  height  of 
230  pm.  As  the  flow  rate  increases,  the 
shear  stress  at  the  location  of  the  sensor 
element  increases  in  a  monotonic  fashion  as 
shown  in  Fig.  15. 

As  the  gas  flows  along  the  flow 
channel,  its  static  pressure  decreases.  Even¬ 
tually  a  limiting  condition  is  reached  where 
the  length  of  the  pipe  cannot  be  increased 
without  altering  the  upstream  conditions; 
that  is,  the  flow  has  choked.  At  that  point, 
dx/dp  is  equal  to  zero.  The  limiting  or  chok¬ 


ing  Mach  number  for  steady  one¬ 
dimensional  isothermal  flow  with  friction  in 
a  constant-  area  duct  is  M=l/Vy.  Thus,  for 
y=1.4,  the  theoretical  choking  Mach  number 
(at  the  channel  exit)  is  0.845.  At  choking  the 
shear  stress  at  the  sensor  location  is  80.55 
Pa,  and  the  Mach  number  at  the  sensor  loca¬ 
tion  is  0.40. 

Fig.  14  shows  the  Mach  number  distri¬ 
bution  along  the  flow  channel  at  the  channel 
height  of  230  pm.  The  distance  from  the 
point  where  the  sensor  is  located  inside  the 
channel  bed  to  the  exit  of  the  channel  is  16.9 
mm.  At  limiting  conditions  for  a  steady  one¬ 
dimensional  isothermal  flow  with  friction  in 
a  constant-area  duct,  the  equation  relating 
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Fig.  14  Mach  number  at  sensor  location  as  a  function 
of  flow  rate  inside  the  calibration  flow  channel. 
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Fig.  15  Shear  stress  variation  with  increase  in  flow 
rate  in  the  calibration  channel. 


the  friction  term  and  the  inlet 

Mach  number  M  is, 


Fig.  16  Typical  response  of  CWRU/Analog  Devices 
MEMS  shear  stress  sensor  at  channel  height  of  230 
pm  (sensor-1). 


sensor  location  implies  a  value 


for  choking.  For  M  =  0.40, 


4-/^nm 

D 


1-yA/2 

yM2 


+  ln  (yM2) 


where,  /  is  the  friction  factor 


1.9682.  At  choking,  the  shear  stress  at  that 
sensor  location  is  80.5  Pa.  The  friction  fac¬ 
tor  f  is  0.0124.  For  Lmax  =  16.9  mm  and  D  = 


/  = 


T 

l/2pV2 


y  is  the  ratio  of  specific 


heats  of  compressed  nitrogen  at  23°C,  D  is 
the  hydraulic  diameter  of  the  channel  [D  =  2 
x  channel  height]  and  Lmax  is  the  length  to 
choking  in  the  channel  from  the  station 
where  the  local  Mach  number  is  M.  Any  in¬ 
crease  in  the  length  would  then  alter  the  up¬ 
stream  conditions.  The  Mach  number  at  the 


0,460  mm,  =  1-84.  This  is  in  close 

agreement  to  the  theoretical  value  of  1.9682. 

Vll-B.  MEMS  Shear  Stress  Sensor 
Calibration  Results 

A  typical  MEMS  shear  stress  sensor  fab¬ 
ricated  at  Advanced  MicroMachines  for  In¬ 
dustry  (AMMI),  now  a  part  of  BF  Goodrich, 
and  packaged  at  Case  Western  Reserve  Uni¬ 
versity  is  mounted  inside  the  flow  channel 
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for  calibration.  The  shear  stress  values  at 
different  flow  rates  for  a  fixed  channel 
height  are  obtained  and  shear  stress  data 
versus  output  voltage  is  plotted.  Figs.  16  and 
17  show  the  measured  output  of  the  mul¬ 
timeter  versus  the  shear  stress  for  two  dif¬ 
ferent  sensors  of  similar  design  at  the  chan¬ 
nel  heights  of  230  pm  and  265  pm  respec¬ 
tively. 

In  Fig.  16,  an  irregularity  in  the  curve  is 
observed  for  a  small  portion  of  the  plotted 
data  before  the  response  behaves  in  a  linear 
fashion.  This  is  due  to  the  stiction  effect. 
The  plot  shows  two  regimes  with  a  different 
sensitivity  where  the  sensor  behaves  line¬ 
arly.  The  sensitivity  of  the  sensor,  at  230  pm 
channel  height,  in  the  regime  I  (up  to  ~  25 
Pa)  is  1.604  mV/Pa  and  in  the  regime  II  (be¬ 
yond  25  Pa)  is  0.57  mV/Pa.  The  linear  re¬ 
sponse  for  regime  I  is  in  the  0.03  V  ~  0.06  V 
range  and  in  the  0.065  V  -  0.1  V  range  for 
regime  n.  The  calibration  accuracy  is  ±2% 
(see  Appendix  A). 


Similar  response  in  the  output  is  ob¬ 
served  in  a  second  sensor  calibration  for  the 
channel  height  of  265  pm.  Fig.  17  is  the  plot 
of  shear  stress  versus  output  voltage  for  the 
calibrated  sensor  a  channel  height  of  265 
pm.  In  regime  I  (up  to  ~  25  Pa),  the  shear 
stress  sensor  exhibits  a  sensitivity  of  about 
0.08  mV/Pa  in  the  linear  operating  range  of 
0.5  mV  to  2.5  mV,  while  in  the  regime  II 
(beyond  25  Pa),  the  sensitivity  of  0.023 
mV/Pa  in  the  linear  operating  range  of  2.5 
mV  to  4  mV  is  observed.  The  calibration 
accuracy  again  is  within  2%. 
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Fig.  17  Typical  response  of  CWRU/Analog  Devices 
MEMS  shear  stress  sensor  for  channel  height  of 
265pm  (sensor-2). 
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It  is  interesting  to  observe  the  similarity 
in  the  transition  from  regime  I  to  regime  II 
in  the  output  of  both  the  sensors  that  were 
calibrated.  The  change  in  the  sensors  sensi¬ 
tivity  starts  when  the  shear  stress  value 
reaches  26  Pa  at  the  sensor  location  and  con¬ 
tinues  till  the  shear  stress  reaches  32  Pa  at 
the  sensor  location.  The  sensitivities  of  both 
the  sensors  at  different  channel  heights  loca¬ 
tions  remains  constant  below  the  shear  stress 
value  of  26  Pa  (regime  I)  and  also  above  32 
Pa  (regime  II). 

The  reason  behind  this  change  in 
sensitivity  is  attributed  to  the  actual  floating 
element  design  on  the  sensor  chip.  In  regime 
I,  the  sensor  is  operating  completely  in  a 
forced-balanced  mode.  The  sensor  element 
is  in  its  null  position  and  the  supporting 
beams  are  not  deflected.  In  regime  n,  the 
floating  element  begins  to  displace  so  that 
part  of  the  restoring  force  is  bome  by  the 
beams  while  the  rest  by  the  capacitive  re¬ 
storing  force.  This  explains  the  reduced  sen¬ 


sitivity  in  regime  II.  Since  both  the  beam 
deflection  and  capacitive  restoration  are  lin¬ 
ear,  the  response  in  regime  II  remains  linear 
but  with  reduced  sensitivity. 

There  is  a  noticeable  jump  in  the  output 
of  Fig.  16  from  0.005  V  to  0.035  V.  This 
jump  is  the  effect  of  stiction,  a  typical  prob¬ 
lem  frequently  encountered  in  the  microfab¬ 
rication  of  microelectromechanical  struc¬ 
tures. 

The  microstructure,  floating  element  in 
the  present  case,  will  be  stuck  to  the  sub¬ 
strate  or  electrostatically  shield  immediately 
beneath  the  surface.  At  very  low  flow  rates, 
the  shear  force  produced  by  the  fluid  flow  is 
not  sufficient  to  overcome  the  stictionai 
forces  beneath  the  surface.  With  the  increase 
in  the  flow  rate,  a  finite  shear  force  that  is 
required  to  overcome  this  stictionai  force  is 
generated  and  the  sensor  element,  which  is 
the  floating  element  in  this  case,  frees  itself 
from  the  substrate.  Once  this  has  occurred 
and  the  shear  forces  dominate  the  stiction 
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forces,  the  floating  element  deflects  in  a  lin¬ 
ear  manner. 

This  problem  could  be  rectified  with  the 
help  of  a  microcontroller  integrated  on  the 
same  chip.  Erasable  Programmable  Read 
Only  Memory  ‘EPROM’  stores  the  offset 
coefficient  and  the  microcontroller  makes 
the  accurate  corrections  so  that  all  the  sen¬ 
sors  with  the  same  specification  exhibit 
similar  output  properties  when  subject  to 
similar  input  conditions.  The  sensors  at  Ana¬ 
log  Devices  Inc.  undergo  a  process  called 
‘Laser  Trimming  of  Resonant  Structures’ 
which  normalizes  the  output  of  all  the  sen¬ 
sor  to  a  standard  output.  This  takes  cares  of 
the  varying  offset  value  in  different  sensors 
of  the  same  kind.  With  the  design  package, 
some  faltering  is  caused  in  the  output  of  the 
sensor  occasionally  due  to  the  imperfections 
and  the  breaks  in  the  gold  line  that  is 
brought  out  of  the  channel  from  the  sensor 
pads  for  external  connections.  This  problem 
can  be  rectified  by  an  enhanced  masking 


process  or  by  thickening  the  gold  line  to 
8000  A°.  The  present  thickness  of  the  gold 
line  is  5000  A°. 

It  has  also  been  observed  that  the  sensor 
shows  photosensitive  characteristics.  Illumi¬ 
nation  causes  drifts  and  wavering  in  the  out¬ 
put  voltage  on  the  order  of  several  micro¬ 
volts.  Fluctuations  in  the  signal  caused  pri¬ 
marily  by  noise,  rotameter  drift,  illumination 
and  measurement  circuit  was  detected  by 
data  acquisition  instruments  and  read  outs 
from  these  instruments  were  corrected  by 
changing  the  offset  value  before  analyzing 
the  collected  data. 

When  the  flow  is  initiated  at  the  start  of 
the  experiment  by  opening  of  the  rotameter 
valve,  a  transient  increase  of  a  few  hundred 
milli-volts  in  the  output  signal  is  observed.  It 
has  a  characteristic  decay  time  of  approxi¬ 
mately  30-35  seconds.  After  the  decay,  the 
output  stabilizes  itself  to  a  steady  drift, 
which  is  then  considered  as  the  offset  value 
for  the  multimeter  and  is  thereupon  included 
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as  a  correction  term.  This  steady  drift  is  con¬ 
sidered  to  be  the  fixed  error  for  the  particu¬ 
lar  experimental  run  and  is  repeatable.  Re¬ 
peatability  can  be  verified  by  obtaining  the 
same  range  of  the  drift  after  closing  and 
opening  the  valve  of  the  rotameter. 

Saturation  of  the  circuit  is  obtained  at 
higher  shear  stresses  because  the  circuit 
cannot  generate  strong  enough  electrostatic 
forces  to  counter  the  shear  forces  caused  by 
the  fluid  flow.  Saturation  level  is  subject  to 
sensor  sensitivity  variation,  which  is  inevi¬ 
table  even  between  mountings  and  position¬ 
ing  of  the  same  sensor. 

VIII.  Uncertainty  Anal  .  sis 

Uncertainty  analysis  refers  to  a  process 
of  estimating  how  great  an  effect  the  uncer¬ 
tainties  in  the  individual  measurements  have 
on  the  calculated  result.  The  square-root-of- 
the-sum-of-the-squares  method,  as  described 
by  Moffat,  is  used  for  computing  an  absolute 
error  of  the  uncertainties  that  are  incorpo¬ 
rated  with  the  calibration.24  In  this  method, 


the  relative  weight  of  each  measurement  is 
taken  into  consideration  as  a  linear  operator 
based  on  its  power  relationship  to  the  value. 
First,  the  fractional  uncertainties  of  the  indi¬ 
vidual  variables  that  are  used  in  the  calibra¬ 
tion  equation  are  measured.  Then,  the  partial 
derivative  of  the  equation  is  taken  with  re¬ 
spect  to  each  component  and  divided  by  the 
total  to  obtain  an  equation  for  how  the  frac¬ 
tional  uncertainties  compound.  Thus,  the 

fractional  uncertainty  ^T/  ;  which  is  also 

defined  as  the  fractional  uncertainty  of  the 
calculated  effective  stress,  can  be  calculated 
from  (9) 
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Applying  the  square  root  of  the  sum  of  the 
squares  rule  to  (9),  the  fractional  uncertain¬ 
ties  of  the  individual  quantities  can  be  ob¬ 
tained; 
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Therefore  combining  all  individual  uncer¬ 
tainties,  we  get  the  fractional  uncertainty  of 
the  calculated  effective  shear  stress*; 


*  The  terms  inside  the  sets  of  parentheses 
are  the  fractional  uncertainties  of  the  individual 
quantities. 
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Each  term  represents  the  contribution  made 
by  the  uncertainty  in  one  factor  to  the  over¬ 
all  uncertainty  in  the  result. 

Taylor  presents  a  systematic  approach  to 
the  identification  and  estimation  of  measur¬ 
ing  system  errors  based  on  manufacturer’s 
specifications 25  The  instruments  that  are 
used  in  the  calibration  process  are  pressure 
transducer,  rotameter,  network  analyzer  and 
multimeters.  Rotameter,  for  the  measure¬ 
ment  of  Mach  number  at  the  sensor  location; 
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Pressure  transducer,  to  calculate  the  pressure 
gradient  along  the  channel  length,  Network 


Analyzer,  to  measure  the  channel  height  and 
multimeters  are  used  for  digital/analog  data 
acquisition. 

The  transducer  calibration  was  found  to 
be  very  well  within  1%  accuracy  as  speci¬ 
fied  on  the  calibration  sheet.  The  inaccura¬ 
cies  of  the  pressure  sensor  from  the  taps 
contribute  in  a  non-linear  fashion  since  the 
gradient  is  the  derivative  of  a  line  based  on  a 
measured  distance.  If  the  pressures  for  the 
curve  fit  are  modified  by  +1%  at  the  maxi¬ 
mum  and  -1%  at  the  minimum  value,  corre¬ 
sponding  to  the  modified  deviance  in  the 
pressure  gradient,  the  pressure  at  the  loca¬ 
tion  of  the  sensor  changes  less  than  0.2%. 


Thus,  maximum  calculated 


dp 

dx 


=  0.002. 


Due  to  the  tendency  of  the  rotameter  to 
drift  during  a  constant  flow,  the  volumetric 
flow  rate  varies  and  leads  to  uncertainty. 
This  would  lead  to  a  slightly  different  pres¬ 


sure  reading  due  to  the  varying  flowrates 
and  the  sensor  output  would  not  correlate  to 
the  recorded  flowrate  and  pressure  gradient. 
To  avoid  this  uncertainty,  a  different  ro¬ 
tameter  with  a  higher  range  from  the  origi¬ 
nal  one  was  used  and  run  at  the  same 
flowrates  for  the  same  test  conditions.  The 
pressure  gradient  and  the  sensor  output  are 
comparable  to  the  one  obtained  from  the 
original  rotameter.  Uncertainty  contributed 
from  this  drift  is  typically  on  the  order  of  5 
to  10%. 

The  relationship  of  the  uncertainty  in  the 
Mach  number  factor  to  that  of  the  total  com¬ 
pressibility  correction  is  obtained  by  taking 
the  derivative  of  the  correction  factor.  The 
compressibility  correction  is  less  than  0.1% 
at  low  shear  stresses  yielding  the  compressi¬ 
bility  correction  fractional  uncertainty  less 
than  0.1%.  The  maximum  compressibility 
correction  fractional  uncertainty  (~  0.5%) 
occurs  at  higher  shear  stresses  where 
M=0.40.  If  the  inaccuracies  in  the  channel 
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height  measurement  and  the  pressure  gradi¬ 
ent  measurements  are  reduced  to  the  level 
comparable  to  that  needed  for  1%  calibra¬ 
tion  accuracy,  the  uncertainty  in  the  com¬ 
pressibility  correction  drops  due  to  the  de¬ 
pendence  on  channel  height  uncertainty  in 
its  calculation.  Thus,  the  uncertainty  in  the 
compressibility  correction  factor  can  be  ne¬ 
glected  for  any  flowrate,  any  shim,  and  any 
reasonable  accuracy  limits  of  calibration. 
The  fractional  uncertainties  of  the  underflow 
shear  stresses  are  neglected  following  the 
same  line  of  reasoning.  The  maximum  un¬ 
certainty  in  the  measurement  of  the  Mach 


number 


(  AM  ^ 


M 


is  0.025. 


Measurements  done  using  the  new  tech¬ 
nique  indicate  a  sensitivity  of  about  ±  3pm, 
which  is  about  ±1.5%  for  a  channel  height 
of  200  pm.  Also,  the  maximum  uncertainty 
in  the  measurement  of  ‘g’,  the  gap  spacing 
between  the  floating  element  and  the  fixed 


substrate 


v  ^  J 


and  ‘t’,  the  thickness  of  the 


floating  element 


v  '  / 


is  0.1.  These  had  neg¬ 


ligible  effect  on  overall  error. 


For, 


'a*' 

dx 

dp 

dx 


=  0.002, 

=0.01, 

(At' 

v  *  J 

v  '  > 

=0.01, 


'A 

v  b  j 


=0.015,  and  AM=0.0 


Therefore  from  (13),  the  total  uncer¬ 
tainty  of  the  effective  wall  shear  stress 


^  At  ^ 

V  Teff  J 


=0.02  at  M=0.4  and  0.015  at  M=0. 


Thus,  MEMS  shear  stress  sensors  are  cali¬ 
brated  to  an  uncertainty  of  ±2%.  The  last  3 
terms  obtained  as  the  uncertainties  from  the 
measurement  of  Mach  number,  the  gap  be¬ 
tween  the  floating  element  and  the  fixed 
substrate  and  the  thickness  of  the  floating 
element  are  negligible  and  are  not  consid¬ 
ered  as  a  major  source  of  uncertainty.  The 
measurements  of  the  channel  height  and  the 
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pressure  gradient  are  considered  to  be  the 
largest  sources  of  error  and  are  a  major  con¬ 
tribution  of  uncertainty  towards  the  overall 
uncertainty  of  the  effective  wall  shear  stress. 

IX.  Summary  and  Conclusion 

A  calibration  flow  channel  and  a  calibra¬ 
tion  procedure  based  on  continuum  isother¬ 
mal  constant-area  compressible  channel 
flow  have  been  developed.  The  calibration 
channel  capabilities  have  been  greatly  in¬ 
creased  and  are  now  able  to  sustain  shear 
stresses  to  80  Pascal.  The  response  of  the 
sensor  was  studied  as  a  function  of  flow  rate 
and  the  channel  height.  A  novel  microstrip 
technique  that  was  implemented  for  measur¬ 
ing  channel  height  was  successful  in  reduc- 
ing  the  inaccuracy  from  ±20  p.m  to  ±3  p.m. 
MEMS  wall-shear-stress  sensors  incorporat¬ 
ing  the  floating  element  technique  with  inte¬ 
grated  electronics  have  been  successfully 
calibrated  with  an  accuracy  of  ±2%. 
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